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Abstract 

In the general two Higgs doublet model with flavor changing neutral Higgs 

couplings, neutral Higgs bosons may decay dominantly via tc or ic final states. 

At the linear collider, e"'"e~ — > h^A^ or H^A^ production processes may result 

in hbtc, W~^W~tc or ttcc (or ttcc) final states. The process 77 — > h^, A^ — > tc 

is also promising, but e^e^ {h^, H^)Z^ tcZ^ is relatively suppressed. 

The possibility of observing like sign lepton pairs, usually the hallmark for 

neutral meson mixing, is quite interesting since mesons do not even form. 
PACS numbers: 14.80.Dq, 14.80.Gt, 12.15. Cc, 13.90.+i 
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Like sign dilepton pair production is the hallmark for heavy neutral meson-antimeson 
mixing. The standard model (SM) predicts rather small mixing effects for mesons containing 
M-type quarks. Furthermore, due to its heaviness, the top quark decays before the T° or T° 
mesons could form. Thus, unlike the b quark case, we do not expect same sign dileptons 
from ti pair production. Effects beyond the standard model are not expected to change this, 
since the Tevatron data is in good agreement with t bW decay dominance expected in 
SM. In this note we report the intriguing possibility of producing like sign top quark pairs 
at linear e~^e~ colliders, within the context of a general two Higgs doublet model (2HDM) 
that possesses flavor changing neutral Higgs (FCNH) couplings f^-^. 

Atwood, Reina and Soni have recently studied FCNH loop induced e~^e~ — > 7*, Z* ^ 
tc transitions at linear colliders and find a rather small rate. They also propose to study 
the tree level s-channel FCNH process fi^fi^ neutral scalars tc. Here, we explore 
FCNH coupling effects in Higgs boson production processes at a 500 GeV e~^e~ Next Linear 
Collider (NLC). We find it to be promising, both for single top tc + X, as well as for the 
more intriguing like sign top pair ttcc final states. 

Let us briefly review the model under consideration. With two Higgs doublets $1 and 
$2; in general one has FCNH couplings. Because of stringent bounds from fi ^ decay, 
K^-K^ and B^-B^ mixings, etc., it is customary |^ to strictly enforce the absence of FCNH 
couplings at tree level. This is readily achieved via some discrete symmetry that allows just 
one source of mass for each given fermion charge 0, much like in SM. However, inspired by 
the quark mass and mixing hierarchy pattern 

mi <^ m2 < m3, 

\Vub? < \Vcb? < \Vus? < 1, 

that emerged since the early 1980 's, Cheng and Sher 0] suggested that low energy flavor 
changing neutral currents could be naturally suppressed, without the need to invoke discrete 
symmetries. Let us elaborate on this observation. 

We shall assume CP invariance throughout the paper, leaving out even the possibility 
of spontaneous CP violation [0. Since both $1 and $2 develop real vacuum expectation 



values (v.e.v.), one can redefine the fields and choose one doublet as the "mass giver", i.e. 
{(pl) = (02) = 0, where v ~ 246 GeV. The usual 2HDM parameter tan/? = vi/v2 gets 

rotated away by the freedom to make linear redefinitions. One readily sees that |^ \/2 Re 0? 
has diagonal couplings. It is, however, not a mass eigenstate. For $2 related fields, we have, 

[uL^^UR + dLi^^^di^ %/2Re0° + [-UL&'^UR + dLi^'^^dn) iV2\m(t>l 

- JiV^t^(")MfiV2 02" + Miye^'^)rfRV2 0+ + H.c, (2) 

where are in general not diagonal, but V^*-^-*^ — ^, since the KM matrix ~ 1. 

At first sight, the Yukawa coupling matrices may appear to be completely general. 
However, in some arbitrary basis where (05) = ^i/a/2, (0°) — quark mass matrices 

consists of two parts, m = m^-^^ + m^'^\ To sustain eq. (1), unless fine-tuned cancellations 
are implemented, one would expect that the off diagonal elements of m^^^ and m^'^\ just like 
m itself, should trickle off as one moves off-diagonal. The rotation (linear redefinition) by 
angle (3 = tan~"'^(t'i/f2) to eq. (2) should not change this property. Hence, data (eq. (1)) 
suggest that C, cannot be arbitrary. In this vein, Cheng and Sher proposed the ansatz 

^ij ~ ^miirLj/v. (3) 

The bonus was that FCNH couplings involving lower generation fermions are naturally 
suppressed, without the need to push FCNH Higgs boson masses to way beyond the v.e.v. 
scale |Tl|. Inspecting eq. (1) again, a weaker ansatz is possible 0, 



i^J = 0{V,^V,^)m^/v. (4) 

According to the mass-mixing pattern, the Cheng-Sher ansatz of eq. (3) corresponds to 
^ij = OiVijVjs) m^/v. Note that in both cases, ^ct is the largest possible FCNH coupling, 
and the associated phenomenology is the most interesting P,p!2|. 

The pseudoscalar = -\/2Ini02 and charged scalar = are already physical 
Higgs bosons, but the neutral CP even Higgs bosons and h'^ are mixtures of -\/2Re05 
and y2Re0^. The mixing angle sin a, a physical parameter, is determined by the Higgs 



potential. In the limit of sin a — * 0, often assumed by various authors ||^, v^Rec/.? 
becomes the "standard" Higgs boson with diagonal couplings, while ^ \/2'Re(f)^ has 
Yukawa couplings as in eqs. (3) or (4), but decouples from vector bosons or charged Higgs 
bosons, just like A^. Our convention for and differs from the minimal supersymmetric 
SM (MSSM) [H, where is taken as the lighter CP even neutral scalar. 

Constraints on the general 2HDM has been studied by various authors. For K*^ and 
mixings, one finds a rather weak bound of 171^0 > 80 GeV, with a more stringent 

bound for . These bounds could weaken, for example, if one uses eq. (4) instead of eq. 



(3). For yU — > 67, an interesting two loop effect dominates over one loop diagrams |jT^. From 



Fig. 4 of ref. [|14l, with ~ 175 GeV, one finds a bound of m/jO > 150 GeV. The bound 



for A^ is weaker since it does not couple to vector bosons and unphysical scalar bosons. 
If sin a — i> 0, the bound would also weaken. A third, less direct constraint on FCNH 
Higgs boson masses is from the recent experimental observation of inclusive b —>■ s'j decays. 
Within the the so-called Model II of 2HDM M (automatically realized in MSSM), the CLEO 



Collaboration gives [1^ a bound of mjj+ > 250 GeV. This bound should weaken for our case 
because of the freedom in ^("''^) as compared to Model II. Inclusion of next-to-leading order 



QCD corrections also tends to soften the bound |[T6|. Thus, we take mH+ ^ 150 — 250 GeV 
as a reasonable lower bound, which is rather consistent with the bounds on FCNH neutral 
scalar bosons. The upshot of our discussion on low energy constraints is that, 

V ~ m(FCNH Higgs) > rut (5) 



is not only reasonable, but quite likely. Although t — c + scalar transitions @JT^ are not 
excluded, we are more interested in Higgs bosons decaying into tc 
We will focus mainly on the mass domain of 

200 GeV < niho^ ao < 2mt ~ 350 GeV. (6) 

We plot in Fig. 1 BR{S° tc + ic) vs. sin^a, for S = h, A and H and m^o = 200, 250, 
300 GeV. A° can decay only via tc and // modes, and can be treated as independent of 



sin^ a. Once m^o is suitably above tc threshold, the tc or tc modes dominate. The behavior 
for at sin a = is similar to A^. However, as sin^ a grows, the VV {V = W, Z) 

partial width grows rapidly, and tc branching ratio becomes rather suppressed. The case 
for is the same as /i" under the interchange of sin^ a < — > cos^ a. The proximity of the 
'^fto j^o = 250 and 300 GeV curves is accidental. 

One might think that the most promising channel for studying FCNH Higgs bosons is 
via the associated production process of e~^e~ Z* ^ H^Z^ and h^Z^ (The W^W~ fusion 
process e^e" vv + is subdominant for the range of eq. (6)). This turns out to be 
not the case. We plot in Fig. 2 the cross section times branching ratio for the signature 
e+e- ^ tcZ^ vs. sin^ a, where S = H, h, again for mgo = 200, 250, 300 GeV. 

For this behaves as x(l — x)/{A{l — x) + ax + 6) where x = sin^ a, A is related to 
the H VV rate in SM, and a, 6 are related to the tc and bb rate. For case one just 
interchanges x ^ 1 — x. Clearly, in the limit of sin a = 0, couples only flavor diagonally, 
while has no production cross section. As sin a grows, because of the mismatch in the 
production and decay process, the effective cross section for tc + Z^ associated production 
remains rather small, with the maximum of 0.43 fb at rrih^fj = 237 GeV and sin^ a (1 — sin^ a 
for H^) = 0.129. This is but a fraction of the total e^e~ — > S^Z^ cross section, which would 
not be easy to observe once one folds in various branching ratios for t or Z decay. 

What is more promising is the e+e~ ^ Z* ^ S^A^ associated production process, where 
S = h, H. In the sin a limit, one has e'^e^ h^A^ only, with cross section similar to 
the H^Z^ mode when phase space is similar. Since in this limit, -/-^ VV, the tc mode has 
a good chance to be the dominant final state for both and We immediately see the 
possibility of our purported ttcc or ttcc final states! 

Let us proceed a bit more systematically. Taking my^o > rriho for illustration (since 
A^ —>■ tc is quite likely to be dominant), we allow rriho to be as low as 100 GeV, with m^o 
in the range of eq. (6), but A^ —>■ h^Z^ is kinematically suppressed or forbidden. For 
numerical illustration, we give in Table | the number of h^A^ events at an NLC with 50 
fb-i integrated luminosity, for sin a = 0, = 200, 250 GeV, m^o = 100, 150, 200 GeV, 
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with the condition that m^o — m^o < 100 GeV. We see that, up to phase space, a few xlO^ 
to 10^ raw events are expected. To determine the number of signal events, one could easily 
fold in branching ratios from Fig. 1, and cos^ a or sin^ a factors for the production cross 
section. These should then be compared with potential backgrounds, the chief ones being 
e+e- ^ tt, Z°Z° events, which are of order 3 x 10^, 4 x 10^ 3 x 10^, respectively. 

Assuming ^ tc + tc is predominant, three modes are of interest: /i° — > bb, tc and VV. 

For m/jO < 2Mw, or when sin^ a <^ 0.1 but the tc mode is suppressed or forbidden, 
/i° bb would be the dominant decay mode. We have a{e'^e'~ h'^A^) x BR(/i° — > 
bb)BR{A^ — > tc + ic) ~ 10 — 20 fb, which is quite sizable compared to the tcZ'^ case. 
However, the bbtc final state may take some effort to identify, since there is just one top 
quark. Although kinematic tricks could be played, but faced with backgrounds that are 
orders of magnitude higher, one would need very good 6-tagging efficiency, and would likely 
need to know m^o beforehand. However, the latter may have to be studied at the NLC 
itself, unless the intermediate mass Higgs search program (via if — > 77 detection) at the 
LHC turns out to be very successful. In any case a detailed Monte Carlo study would be 
necessary to determine whether this mode can be fruitfully studied. 

For sin^a <S 0.1, as can be seen from Fig. 1, /i — > is hkely dominant in the mass 
range of eq. (6). We find a{e+e- h^A°) x BR{h° tc + ic)BR{A° tc + ic) ~ 10 fb, 
which is shghtly smaller than the previous case because of phase space. Clearly, 50% of this 
cross section goes into ttcc or ttcc final states, which is again larger than the tcZ° case. We 
expect typically of order 250 such events. Folding in the semileptonic branching ratio, one 
expects ~ 12 events in the signal of 

e+e- — > + + 4j, (7) 

where the 4 jets have fiavor bbcc or bbcc. Thanks to the large top quark mass, this distinctive 
signature has seemingly no background. In contrast, if one allows only one top to decay 
semilcptonically, or if one tries to probe the equivalent number of ttcc events, the single i + 
u+Gj or opposite sign dilepton I'^i'^ + uu+Aj signatures would be swamped by tt or W~^W~ 
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production background, which are orders of magnitude higher. In particular, standard 
e^e^ ti pair production with hard gluon radiation may be especially irremovable. Since 
the effect demands m/jO + m^io > 400 GeV, one is phase space limited at a 500 GeV NLC If 
the center of mass energy of the NLC could be increased to 600 GeV or so, possible phase 
space suppressions for producing e^e'^ — » h^A^ — > ttcc + ttcc could be relieved. 

For sin^ a ^ 0.1 and rriho > 2Mw, W^W~ and Z^Z^ decays are (pre) dominant, 

with tc no more than 10%, and like sign top quark pair final states become no longer 

visible. The case remains interesting, however, since W^W'tc or Z^ZHc final states at 
the 5 — 10 fb level (taking into account both and H^) are still quite conspicuous. One 
effectively has W^W~W^hc in final state, which again has little background. In particular, 
one could still have like sign dilepton pairs as in eq. (7). With m/jo or m^/o known from the 
LHC, this decay mode could be studied at the NLC in complete detail. 

It is also possible to produce Higgs boson via the 77 — > process 0. We note that 
only couples to fermions, hence its effective coupling to photons is smaller than the SM Higgs 
boson. For if sin a is very small, the case is again similar. As sin a grows, the effective 
coupling would quickly become dominated by vector bosons and the production cross section 
could be larger. This is, however, offset by the reduction in tc branching ratio. Since 

77 ^ 5'° — >• tc should have little background {e.g. tc — ^ W^hc can be distinguished from 
W^W~ via 6-tagging since BR(iy~ hq) < 10^'^), the number of events expected is 

Ar(77 S° ^ tc + ic) = 47r2r(5° 77) BR(5° ^tc + ic)/m% 

X (1 + (AA')) {dC^^/dm^^)\m-,^=mgo- (8) 



It is possible [|T^ to tune photon polarizations to have (AA') ~ +1 and effective luminosities 
close to the the e^e~ mode [i.e. ~ 50 fb^^). If such is the case, then one expects 10^ — 10'^ raw 
events, which should make tc detection possible if the branching ratio is not too suppressed. 
Note that the corresponding number of 77 — > W^W~ pairs is at the 10'^ — 10^ level. 

We now compare our results with that of Atwood, Reina and Soni. For e"'"e~ — > 7*, Z* 
tc via FCNH loop effects, they find i?*'' = cr(e+e" ^ tc + tc)/(T(e+e~ — * 7* A*'*'/^") ~ 



few X 10~^ which amounts to less than 0.1 event for a 500 GeV NLC with 50 fb~^ 
integrated luminosity. In case and are heavier than the range of eq. (6), the loop 
induced cross section also goes down by another order of magnitude 0. Thus, this process 
is unlikely to be observable at the NLC. Although phase space favored, loop suppression in 
this case is too severe. For yU+yU^ — > h^, tc + ic, the process occurs at tree level and 

has a sizable cross section [0. But in the limit of sin a — 0, would also not decay via the 
VV mode, just like A^. A rather fine-stepped energy scan would then be needed because of 
the narrowness of the and A^ width. Together with the technological uncertainty for a 



high energy, high luminosity fi^fj^^ collider |[T9|, this process might be less straightforward 
to study than at the NLC, including the 77 collider option via 77 h^, A^. 

As stated in the Introduction, the signature of like sign top pair production is rather 
analogous to observing T^T^ or T^T^ pairs via T^-T^ mixing. With top mesons not even 
forming, however, it is the associated production of h^A^ pairs, which each subsequently 
decay equally into tc or ic final states, that circumvents the usual condition of associated 
production of ti (or, W^W~) pairs in most processes. Since and A^ contribute to B-B 
mixing, in a sense the like sign top pair production effect is related to neutral meson mixing 
phenomena. We know of no other way to make tt or tt pairs in an e^e~ collider environment. 

In summary, within a general two Higgs doublet model with FCNH couplings, , A^ ^ 
tc+ic could be the dominant decay mode. The most intriguing consequence is the possibility 
of detecting like sign top pair production via e^e~ h^A^ ttcc or ttcc, while single top 
tcbb or tcW^W~ production are also detectable. In contrast, the e"'"e^ Z^) — > tcZ^ 

process is rather suppressed and not competitive. The number of events, hence the FCNH 
Higgs boson mass reach, could be extended if the collider energy is higher. It is also possible 
to study single FCNH Higgs production via 77 A° ^ tc + tc. Since the neutral 

Higgs bosons of minimal supersymmetric standard model couple to fermions in a flavor 
diagonal way, the observation of FCNH signals would rule out MSSM. We urge experimental 
colleagues to study signal vs. background issues carefully. 
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FIGURES 

FIG. 1. BR{S^ tc + ic) vs. sin^a for S = h (dash), A (dotdash) and H (solid) and 
m-so = 200, 250, 300 GeV (from bottom to top). 

FIG. 2. c7(e+e- ^ S^Z°)x BR(5° ^ tc) ?;s. sin^ a, ioi S = h (dash), (sohd) and mgo = 200 
(lower curve), 250 (upper curve), 300 (middle curve) GeV. 
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TABLES 



TABLE I. Number of e+e' /i°A° events for sina = at NLC with JjCdt = 50 fb^^ 



m^o (GeV) 


200 


200 


250 


200 


250 


irifio (GeV) 


100 


150 


150 


200 


200 




1160 


900 


490 


520 


200 
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